The Shear Viscosity and Thermal Conductivity of Nuclear Pasta 
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We calculate the shear viscosity 77 and thermal conductivity k of a nuclear pasta phase in neutron 
star crusts. This involves complex non-spherical shapes. We use semiclassical molecular dynamics 
simulations involving 40,000 to 100,000 nucleons. The viscosity r/ can be simply expressed in terms 
of the height Z* and width Aq of the peak in the static structure factor Sp{q). We find that 
r] increases somewhat, compared to a lower density phase involving spherical nuclei, because Z* 
decreases from form factor and ion screening effects. However, we do not find a dramatic increase 
in ri from non-spherical shapes, as may occur in conventional complex fluids. 

PACS numbers: 26.60.-c Nuclear matter aspects of neutron stars, 97.60.Jd Neutron Stars, 66.20.Cy Theory 
and modeling of viscosity of liquids..., including computer simulations 



I. INTRODUCTION 

The viscosity of dense matter can dampen the oscilla- 
tions of neutron stars. High multipolarity p-mode oscilla- 
tions may impact the pulse shapes of some radio pulsars 
[l|]. For p-modes, the primary restoring force is the pres- 
sure, and the modes may be damped by the shear viscos- 
ity of the neutron star crust [2]. Inertial or r-modes of 
a rotating neutron star may radiate gravitational waves 
and could limit the spin period iS]. For these modes the 
coriolis force provides the primary restoring force and 
they may be damped by the shear, and or bulk, viscosi- 
ties of dense matter through out the star. 

Flowers and Itoh [4| and Nandkumar and Pethick [5j 
have calculated the shear viscosity of neutron star crust 
matter in the liquid phase. While recently, Chugunov 
and Yakovlev have calculated the shear viscosity of both 
the solid and liquid phases [2[- They assume the viscos- 
ity is dominated by the momentum carried by degenerate 
electrons and that this is limited by electron-ion scatter- 
ing. At high densities in the neutron star core, the shear 
viscosity comes from momentum carried by nucleons and 
this is limited by NN scattering. See for example Ref. 
i. 

The thermal conductivity of the neutron star crust 
helps determine the temperature profile of the star and 
is very important for the cooling time of the crust. Re- 
cently, crust cooling times have been observed for two 
neutron stars after extended outbursts Q- We have cal- 
culated the thermal conductivity of the outer crust in- 
cluding the possible role of impurities [g]. The thermal 
conductivity could also be important for setting the ig- 
nition conditions for carbon superbursts [^|. 

At intermediate densities, around lO^"* g/cm'^ just be- 
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low nuclear density, matter may form complex nuclear 
pasta phases [lO|. Competition between short range at- 
tractive nuclear, and long range repulsive Coulomb, in- 
teractions can lead to clusters with many different non- 
spherical shapes including long rods or flat plates. Be- 
cause pasta may form at high densities, it could represent 
as much as half the mass of the neutron star crust. The 
complex shapes in nuclear pasta have sizes of tens of Fer- 
mis. This is comparable to the wavelength of neutrinos 
in Supernova explosions. Therefore, coherent scattering 
from the nuclear pasta shapes may impact the opacity of 
Supernova neutrinos |ll| . 

In the laboratory, complex fluids with large non- 
spherical molecules can have shear viscosities that are 
much larger than for simple fluids. How might the com- 
plex pasta shapes influence the shear viscosity of nuclear 
pasta? Is there a drastic increase in viscosity that is sim- 
ilar to that seen for laboratory complex fluids? We are 
not aware of any previous calculations of the shear viscos- 
ity of nuclear pasta. Although Chugunov and Yakovlev 
have calculated the shear viscosity for similar densities 
[2I, they apparently assumed spherical nuclei. 

In this paper, we calculate the shear viscosity and ther- 
mal conductivity using a semiclassical model for the nu- 
clear pasta [I4I and molecular dynamics (MD) simula- 
tions with from 40,000 to 100,000 nucleons [l3|. Our 
goal is to gain some simple insight into how transport 
properties, such as the shear viscosity or thermal con- 
ductivity, depend on the sizes and shapes of the clusters. 
Furthermore, we wish to know how transport properties 
might change as one goes from isolated spherical nuclei 
at low densities, to complex pasta phases, and then to 
uniform nuclear matter at high densities. One important 
advantage of our semiclassical model is that it can be 
directly applied in all of these density regimes. 

We wish to emphasize the importance of simple qual- 
itative insight. Astrophysicists request very detailed 
transport properties, such as the shear viscosity, even 
though there are large remaining uncertainties in much 



more basic dense matter properties. It may be unre- 
alistic to predict the exact sizes and shapes of the pasta 
clusters. Instead, we hope to gain qualitative insight into 
how pasta shapes may impact the shear viscosity or other 
properties. 

The manuscript is organized as follows. In Sec. UIl we 
review the simple semiclassical model, explain the cal- 
culation of the shear viscosity, and provide some details 
of the molecular dynamics simulations. In Sec. IIIII we 
present results for the static structure factors and use 
these to calculate the shear viscosity and thermal con- 
ductivity. Finally, conclusions and future directions are 
presented in Sec. HVl 



of a sum of two-body interactions 



II. FORMALISM 



In Subsection III Al we review our semiclassical model 
that can describe nuclear pasta phases in a flexible way. 
Next, in Subsection lll BI we describe the calculation of the 
shear viscosity based on molecular dynamics simulations. 



A. Semiclassical Nuclear Pasta Model 

We briefly describe our semiclassical model that while 
simple, contains the essential physics of competing in- 
teractions consisting of a short-range nuclear attraction 
and a long-range Coulomb repulsion. This competition 
can generate complex pasta shapes. The impossibility 
to simultaneously minimize all elementary interactions is 
known in condensed-matter circles as frustration. The 
complex physics of frustration, along with m any other 
details of the model, may be found in Refs. [i3j[3- Here 
only a brief review of the most essential features of the 
model is presented. We model a charge-neutral system 
of electrons, protons, and neutrons. The electrons are 
assumed to be noninteracting and thus are described as 
a degenerate free Fermi gas at a number density iden- 
tical to that of the protons {i.e., Ue = Pp). The nucle- 
ons, on the other hand, interact classically via a nuclear- 
plus-Coulomb potential. However, the use of an effective 
temperature and effective interactions are used to simu- 
late effects associated with quantum zero-point motion. 
While simple, the model displays the essential physics of 
frustration, namely, nucleons clustering into pasta but 
the size of the clusters limited by the Coulomb repulsion, 
in a transparent form. Moreover, one may study the evo- 
lution of the system through the low density, pasta, and 
high density phases within a single microscopic model. 
Finally, the model facilitates simulations with a large 
numbers of particles, a feature that is essential to esti- 
mate and control finite-size effects and to reliably study 
the long wavelength response of the system. 

The total potential energy V^ot of the system consists 
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where the "elementary" two-body interaction is given as 
follows: 

V{i,]) = ae~^"^'^ +[h + cT,{i)T,X3)\e-'-'^^'^'' + VS,j) ■ 

(2) 
Here the distance between the particles is denoted by 
rij = \Yi — Tj\ and t^ represents the nucleon isospin 
projection (r^ — 4-1 for protons and r^ = —1 for neu- 
trons). The two-body interaction contains the charac- 
teristic intermediate-range attraction and short-range re- 
pulsion of the nucleon-nucleon force. Further, an isospin 
dependence has been incorporated in the potential to en- 
sure that while pure neutron matter is unbound, sym- 
metric nuclear matter is appropriately bound. Indeed, 
the four model parameters (a, 6, c, and A) introduced in 
Eq. ^ have been adjusted in Ref. [l2| to reproduce the 
following bulk properties: a) the saturation density and 
binding energy per nucleon of symmetric nuclear matter, 
b) (a reasonable value for) the binding energy per nucleon 
of neutron matter at saturation density, and c) (approx- 
imate values for the) binding energy of a few selected 
finite nuclei. All these properties were computed via a 
classical Monte Carlo simulation with the temperature 
arbitrarily fixed at 1 MeV. The parameter set employed 
in all previous and present calculations is displayed in 
Table H] Finally — and critical for pasta formation — a 
screened Coulomb interaction of the following form is in- 
cluded: 



VS,j)^—e-''^''^rp{i)Tp{3), 
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where Tp = {1-\-Tz)/2 and A is the screening length that 
results from the slight polarization of the electron gas. 
The relativistic Thomas-Fermi screening length is given 

by 



A = ^ ( k-pxlkl+ml 
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(4) 



where nie is the electron mass, the electron Fermi mo- 
mentum has been defined by fcp = (STr^rie)^''^, and n^ is 
the electron density [l^l . Unfortunately, while the screen- 
ing length A, defined above, is smaller than the length 
L of our simulation box, it is not significantly smaller. 
Hence, following a prescription introduced in Ref. [13| 
in an effort to control finite-size effects, the value of the 
screening length is arbitrarily decreased to A = 10 fni. 
This decrease in the screening length could slightly in- 
crease the size of clusters because it somewhat reduces 
Coulomb repulsion. However, we do not expect this to 
be a large change because A is still larger than, or com- 
parable to, the size of clusters. In the future, MD simula- 
tions should be performed using the full physical screen- 
ing length. 



The simulations are carried out with both a fixed num- 
ber of particles A and a fixed density p. The simula- 
tion volume is then simply given by T^ = A/ p. To mini- 
mize finite-size effects periodic boundary conditions are 
used. To carry out molecular dynamics simulations the 
trajectories of all of the particles in the system are de- 
termined by simply integrating Newton's laws of motion, 
albeit for a large number of particles (up to 100,000 in 
the present case) using the velocity- Verlet algorithm [lj| . 
To start the simulations, initial positions and velocities 
must be specified for all the particles in the system. The 
initial positions are randomly and uniformly distributed 
throughout the simulation volume while the initial ve- 
locities are distributed according to a Boltzmann distri- 
bution at temperature T. As the velocity- Verlet is an 
energy — not temperature — conserving algorithm, kinetic 
and potential energy continuously transform into each 
other. To prevent these temperature fiuctuations, the 
velocities of all the particles are periodically rescaled to 
ensure that the average kinetic energy per particle re- 
mains fixed (3/2)fcBT. 

In summary, a classical system has been constructed 
with a total potential energy given as a sum of two- 
body, momentum-independent interactions as indicated 
in Eq. ([2|). Expectation values of any observable of in- 
terest may be calculated as a suitable time average using 
particle trajectories generated from molecular dynamics 
simulations. 

We comment on the classical nature of our simulations. 
Correlations from Fermi statistics are not explicitly in- 
cluded. However, some effects of Pauli exclusion are im- 
plicitly included by incorporating short range repulsion in 
Eq. ^ and adjusting the parameters in Table[T]to repro- 
duce the saturation density and binding energy of nuclear 
matter. In this paper, we focus on the electron-proton re- 
sponse at long wavelengths. This is dominated by corre- 
lations between clusters. These clusters are heavy, since 
they involve many nucleons, and their thermal de Broglie 
wavelengths are much shorter than the inter-cluster spac- 
ing. Therefore, we expect our classical approximation to 
accurately reproduce the long wavelength response. 



B. Shear Viscosity 

Following Chugunov and Yakovlev Q, we assume the 
shear viscosity rj is dominated by momentum carried by 
electrons and that this is limited by electron-ion scatter- 
ing, 

V = 7^77^ -JT- ■ (5) 

Chugunov et al work in ion coordinates with Zion the 
ion charge, the electron Fermi velocity is uf ~ 1, and 
Tie is the electron density. The Coulomb logarithm Aei 
describes electron ion scattering , 

'■2'=^ dq F{qf 
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Here e^q) is the static longitudinal dielectric function of 
the electron gas [lal (^(q) ~ 1 except at low momentum 
transfer q). The static structure factor Sion{q) describes 
correlations between ions and can be calculated from the 
ion density-density correlation function, see below. The 
nuclear form factor F{q) is. 



nq) 
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dPre'"^" Pp{r) , 
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where Pp{r) is the proton density inside one ion {F{q = 
0) = 1). 

As the density increases, and the nuclei start to 
strongly interact to form complex pasta shapes, one 
may no longer be able to effectively use ion coordinates. 
Therefore we work directly in the nuclcon coordinates 
instead of using the ion coordinates. At low densities 
the nucleons form clusters in our model that are equiva- 
lent to nuclei. As a result our calculation in the nucleon 
coordinates, at low densities, will reproduce a calcula- 
tion in the ion coordinates. We replace the Coulomb 
logarithm describing electron-ion scattering with one de- 
scribing electron- proton scattering A^p, 
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The electron-proton Coulomb logarithm is, 
'^"^ dq f.^iq)^ 
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Here fsn (9) is the single nucleon form factor (charge dis- 
tribution of the proton) that we approximate fsn{q) — 1- 
The static structure factor Sp{q) describes correlations 
between protons and is calculated from the density- 
density correlation function 



Sp{q) = {ppiq)* Pp{q)) 

Here Pp{q) is the proton density 

Pp{q) 
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for a simulation with Np protons at positions '^iii) and 
the statistical average in Eq. [Tl] is calculated as an aver- 
age over the time t in the MD simulation. 

Our model automatically includes the correlations be- 
tween nucleons that go into forming ions of charge Zion- 
Furthermore, the model also includes the distributions 
of protons inside a nucleus that is described by the form 
factor F{q). In the limit where Zjo„ protons correlate 
into each cluster and there are only weak correlations 
between the clusters one would have 



Sp{q) w ZionFiqY 



(13) 



TABLE I: 


Model parameters used in the calculations. 


a 


h c k 


110 MeV 


-26 MeV 24 MeV 1.25 fm^ 



TABLE 11: Simulation parameters and shear viscosity results: 
the baryon density is p, the total number of particles in the 
simulation A, and Ttot is the total simulation time. The height 
of the peak in Sp{q) is Z* , at location q* and full width Ag. 
The Coulomb logarithm, Eq. 1101 is hep and Z* ^q/q* provides 
a simple approximation to Aep, see text. The shear viscosity 
(in units of 10^" Pascal-seconds) is r). 



A 



fm" 



Ttot 
fm/c 



Z* q* Aq 
fm^^ fm^^ 



Z'Aq 
9* 



Ae 



lO^'^ Pa-s 



0.01 40000 1.3 X 10" 
0.025 100000 52000 
0.05 100000 28000 



34.4 0.28 0.108 13.3 13.7 5.3 
38.6 0.31 0.077 9.5 10.4 17 
20.3 0.38 0.090 4.9 5.8 63 



MD Simulations 



The MD simulations of ref. [l^ were for a temperature 
of 1 MeV, a proton fraction of Yp = 0.2, and baryon den- 
sities p of 0.01, 0.025, and 0.05 fm'^. See TableHH Note 
that the proton fraction Yp = 0.2 is intermediate between 
the higher proton fraction expected in a supernova core 
during the in fall phase and the lower proton fraction ex- 
pected for neutron star crust in beta equilibrium. Figure 
[T] shows the 0.03 fm~'^ iso-surface of the proton density 
for one configuration of the p = 0.01 fm^"^ simulation. At 
this density, all of the protons and most of the neutrons 
are clustered into nuclei. There is also a low density neu- 
tron gas between the clusters which is not shown. Figure 
[5] shows the proton density at the higher baryon density 
p = 0.025 fm'^. Now the clusters are larger (larger mass 
nuclei) and are closer together. The nuclei are larger be- 
cause the higher density electron gas cancels more of the 
Coulomb repulsion. This allows nuclei to form with more 
protons. 



and 



so that 



Sioniq) ~ 1 



ZionSion{q)F{qY ~ Sp{q) . 



(14) 
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In this paper we are interested in transport properties 
of electrons. These depend on the proton static structure 
factor Sp{q). However for comparison, we also consider 
the static structure factor for neutrons S'„ (q) . This is im- 
portant for neutrino transport. Neutrinos in a supernova 
scatter from the weak charge density. In the Standard 
model the weak charge of a neutron is much larger than 
that of a proton. Therefore neutrino scattering can be 
described with the static structure factor for neutrons 
Sniq), see for example [l3||ll|. 



Sn{q) = {pn{q)* Pn{q)) , 

Here Pn{q) is the neutron density 



(16) 




FIG. 1: (Color online) The 0.03 fm~ proton density isosur- 
face for one configuration of 40,000 nucleons at a density of 
0.01 fm"''. The simulation volume is a cube 159 fm on a side. 



N„ 



Pn{q) 






j^-<^i{t) 



(17) 



where the sum runs over Nn neutrons at positions ri{t). 
We will calculate Sp{q) and Sn{q) in Section [1111 



III. RESULTS 

In this section we present MD simulation results for 
static structure factors, shear viscosity, and thermal con- 
ductivity. 



Figure[3]shows the proton density for one configuration 
of 100,000 nucleons at a density of 0.05 fm~'^. The clus- 
ters are now seen to have very elongated shapes. The low 
density neutron gas between these clusters is not shown. 
These elongated spaghetti like shapes are one example 
of a nuclear pasta phase. Note that the shapes are not 
straight rods but instead they bend and branch inside 
the simulation volume. Because of the periodic bound- 
ary conditions, the shapes extend out one side and back 
in another side of the simulation volume. 




FIG. 2: (Color online) The 0.03 fm^ proton density isosur- 
face for one configuration of 100,000 nucleons at a density of 
0.025 fm~'^. The simulation volume is a cube 159 fm on a 
side. 




FIG. 3: (Color online) The 0.03 fm~^ proton density isosur- 
face for one configuration of 100,000 nucleons at a density of 
0.05 fm~^. The simulation volume is a cube 126 fm on a side. 



B. Static Structure Factor Results 

We calculate the static structure factors Sp{q) for pro- 
tons and Sn{q) for neutrons at p = 0.01 fm~^ by av- 
eraging over 16000 configurations of the 40000 nucleons 
where each configuration is separated from the next by a 
simulation time of 20 fm/c. Note that calculating the 
static structure factor by averaging over a finite time 
16000 X 20 = 320000 fm/c may slightly under estimate 
the height of the peaks in Sp{q) and Sn{q) compared to 



averages over much longer times. We consider momen- 
tum transfers q, 



27r 
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with L the length of the simulation volume and Ux, n 
are integers. We average over the directions of q. Figure 
m shows Sp{q) and 5„(g) at p = 0.01 fm~^. 
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FIG. 4: (Color online) The static structure factor Sp{q) for 
protons (solid black line) , and for neutrons Sn (g) (dotted red 
line) versus momentum transfer g at a baryon density of 0.01 
fm-^ 



Our results for Sn{q) agree well with the original re- 
sults for S{q) presented in ref. \\M,- Note that in ref. 
[13| . S[q) was calculated from the Fourier transform of 
the radial distribution function g{r) (because this was 
available). However this involved some uncertainty from 
the assumed form of g{r) at large distances r. Our re- 
sults for Sp{q) are new. This was not calculated in ref. 

The static structure factor Sn(q) is larger than Sp{q) at 
small q. This is because long wavelength fluctuations in 
the neutron density need not feel an electrostatic restor- 
ing force. Therefore they are not strongly screened. In 
contrast all long wavelength fluctuations in the proton 
density are strongly screened and this makes Sp{q) small 
at small q. This has been discussed in ref. [l^. See also 
ref. [13. 

At intermediate g, there is a large peak in both Sp{q) 
and Sn{q)- This corresponds to elastic electron- nucleus 
(for Sp{q)) or neutrino- nucleus (for Sn{q)) scattering. 
The cross section for coherent elastic electron-nucleus 
scattering goes like the atomic number squared Z'^ . How- 
ever, the structure factor Sp{q) is normalized per proton 
(instead of per ion). Therefore, the cross section per pro- 
ton is proportional to Z. One can think of the height of 



the peak in Sp{q) as being the effective number of pro- 
tons Z* that the electron scatters from. This effective 
number of protons Z* is reduced from the actual number 
of protons in a nucleus because of the form factor F{q), 
Eq. [7l The form factor leads to a large reduction at 
high q. In addition Z* is reduced by the screening effects 
of other ions. This reduction is very large at small q. 
Finally correlations between ions can increase Z* some- 
what, see ref. [l^. This corresponds to Sioniq) > 1- We 
characterize the peak in Sp{q) by its height Z*, its full 
width at half maximum Ag, and its location q*. These 
values are illustrated in Fig. [S] and collected in Table HIl 



higher at p = 0.025 fm-^ than at p = 0.01 fm-^. This 
is because screening from the other ions is more effec- 
tive and extends to higher q at this higher density. This 
limits the peak in Sp{q) from the left (for low q). Also 
the nuclear form factor F(q) falls more quickly with q 
because the nuclei are now larger. This limits the peak 
in Sp{q) from the right (for high q). As a result the peak 
is narrower (has a smaller Ag) at a density of 0.025 than 
at 0.01 fm-3. 
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FIG. 5: The static structure factor Sp{q) for protons versus 
moinentuiii transfer g at a baryon density of 0.01 fm""^ (note 
the linear scale). The height of the peak in Sp{q) is Z* , the 
full width at half maximum is Ag, and the peak position is 

q*- 



s2^ 









A 






10 


r 


/ 


'V 




— Protons 
Neutrons 


1 


r 


/ 




^xv 


v_^: 


H 


■/ 










) 


0.2 


0.4 


0.6 


0.8 1 



q (fm" ) 

FIG. 6: (Color online) The static structure factor Sp{q) for 
protons (solid black line) , and for neutrons S„ (g) (dotted red 
line) versus momentum transfer g at a density of 0.025 fm~^. 



The system is neutron rich. Therefore one might ex- 
pect the peak in S'„(g), representing an effective num- 
ber of neutrons N* , to be larger than the peak in Sp{q). 
However this is not the case. A neutrino will scatter from 
neutron density contrasts. The more or less uniform neu- 
tron gas, that is also present outside the nuclei, leads to 
a reduction in contrast. One only scatters from the dif- 
ference in inside to outside neutron density. As a result 
the peak in 5'„(g) is smaller than the peak in Sp{q) even 
though the nuclei are very neutron rich. 

We calculate the static structure factors Sp{q) and 
Sn{q) at p = 0.025 fm~^ by averaging over 10500 config- 
urations of the 100000 nucleons, of the simulation shown 
in Fig. [51 Each configuration is separated from the next 
by 20 fm/c. Figure [5] shows Sp{q) and S'„(g) at p — 0.025 
fm^'^. Again there is a large peak in Sp{q) that is slightly 
higher than the peak at p = 0.01 fm~'^, in Fig. |4l The 
nuclei at p — 0.025 fm~'^ have about twice the number of 
nucleons A w 200 as those at p = 0.01 fm~^ which have 
A w 100, see ref. [ll|. Nevertheless Z* is only slightly 



Finally we calculate Sp{q) and ^^(g) at a density 0.05 
fm^'' by averaging over 19400 configurations of 100000 
nucleons. Again each configuration is separated by 20 
fm/c. At this density the system is in a nuclear pasta 
configuration with complex rod-like shapes, see Fig. [S] 
Figure [7] shows ^^(g) and S'„(g) at 0.05 fm~^. There is 
still a peak in Sp{q). However it is not as high as at 
lower densities. Each piece of pasta contains a very large 
number of nucleons. However, they are packed closely 
together and this allows screening from other pieces of 
pasta to be very effective. This greatly reduces Sp{q). 
Figure [7] shows no qualitative effect from the non-sperical 
pasta shapes. Instead Sp{q) has very similar behavior for 
round nuclei in Figs. I4|6I or long rods in Fig. [7l Note 
that we have averaged over the directions of q. The sys- 
tem can have large pasta shapes with very many nucle- 
ons. However, it is an important result of this paper that 
screening and form factor effects limit the number of pro- 
tons that one can coherently scattered from, for any one 
momentum transfer. 
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FIG. 7: (Color online) The static structure factor S-p{q) for 
protons (solid black line), and for neutrons Sn{q) (dotted red 
line) versus momentum transfer g at a density of 0.05 fm~^. 
At this density the system is in a nuclear pasta phase. 



C. Shear viscosity results 



We use results for Sp{q) in Figs. l4l6l7l to calculate the 
Coulomb logarithm Aep, in Eq. IIOI and the viscosity rj, 
from Eq. [^1 These results are collected in Table HIl To 
gain insight into the results we approximate Sp{q) in Eq. 
[To] with a rectangle of height Z*, width Aq, and centered 
at q* so that Agp ~ AqZ* /q*. This yields a very simple 
approximation for r/, 
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Therefore, -q follows from the parameters Z*, q* and Aq 
that characterize Sp{q). 

The position of the peak q* is expected to increase 
slightly with density as the clusters are packed more 
closely together. Perhaps more interesting is the behav- 
ior of the width Aq with increasing density. The width 
describes the range of momentum transfers over which 
e — p scattering is effective. This is limited for low q by 
ion screening and at high q by the nuclear form factor. 
Screening can be effective for wavelengths longer than 
the distance between ions, while the form factor is small 
for wavelengths much smaller than the size of a nucleus. 
Therefore we have the following very approximate way 
to think about Aq. 



Aq 



1 



Cluster size Distance between clusters 



(20) 



In going from a density of 0.01 to 0.025 fm^'^ the distance 
between clusters (ions) decreases. In addition the cluster 
size increases as the composition changes from nuclei with 



AKilOOto AKi 200. Thus Aq "gets squeezed from both 
sides" and it is smaller at a density of 0.025 than it is at 
0.01 fm-3 

The effective number of protons Z* is close to the 
actual number of protons in a nucleus Z at low densi- 
ties. However as the density increases, screening and 
form factor effects greatly reduce Z* so that Z* ^ Z. It 
appears that simple geometric constraints keep Z* rela- 
tively small even if Z is large. Indeed Table |ll] shows that 
Z* decreases as one goes from p = 0.025 to 0.05 fm~'^ 
and the system changes from isolated nuclei to very long 
pasta shapes that contain large numbers of protons. 

Our results for t] are collected in Table [TTl Note that 
these values use the full integration in Eq. [10] instead of 
the approximate form in Eq. 1191 However the difference 
is small. Chugunov et al. [4I , in calculations for spherical 
nuclei, find the same order of magnitude for 77. However 
they only present results for smaller temperatures and 
proton fractions. This confirms our finding that non- 
spherical pasta shapes do not greatly change the shear 
viscosity. 

There can be large uncertainties in the pasta shapes or 
even the densities where pasta phases are present. How 
might these uncertainties impact 77? The viscosity, for 
whatever complex shapes are present, can be calculated 
in terms of the simple effective parameters Z* , Aq, and 

q ■ 

We comment on our use of the Born approximation 
to describe e-ion scattering. Second or higher order 
Born corrections will modify our results somewhat, see 
for example [l^, but should not change them qualita- 
tively. Note, there is some ambiguity in determining the 
best value of the nuclear chage Z to use in calculating 
Coulomb distortions in the pasta phase. Perhaps one 
could use Z*. 

We now discuss the temperature dependence of rj. We 
have performed calculations at the relatively high tem- 
perature of T = 1 MeV, because this was the tempera- 
ture of our MD simulations in ref. [lj|. In general we 
are interested in rj for lower temperatures. Equation 1191 
shows that rj depends on T only through the tempera- 
ture dependence of the parameters Z*, Aq, and g*. We 
expect the physical size of the pasta clusters to be sim- 
ilar at lower T to their size at T = 1 MeV. However, 
screening will be more efficient at lower T. This could 
decrease Z* and somewhat increase ry. The temperature 
dependence of rj should be studied in future work with 
MD simulations at lower temperatures. 

Finally, we comment on neutrino contributions to rj. 
At high temperatures, for example during a supernova, 
momentum carried by neutrinos should dominate the 
shear viscosity because of their very long mean free path. 
However at low temperatures, neutrinos will no longer re- 
main in equilibrium with matter and their density will be 
very low. As a result, 77 will be dominated by the electron 
contributions. 
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TABLE III: Thermal conductivity k results. The Coulomb 
logarithm, Eq. [551 is A^p. 



P (fP^' 



a;: 



K ( eig/K cm s) 



0.01 


16.1 


1025 


11.7 


0.05 


6.3 



2.5 X 10^° 
4.7 X 10^° 
1.1 X 10^^ 



D. Thermal Conductivity 

The thermal conductivity k of nuclear pasta can be 
calculated in an almost identical way to 77. We assume 
the conductivity is dominated by the energy carried by 
electrons and that this is limited by electron-ion scatter- 
ing. Following Nandkumar and Pethick i5| we write, 



12a2A-p 



(21) 



where the Coulomb logarithm AJL is slightly different 
from the corresponding one A^p for the shear viscosity. 



A*" = 

ep 



2kf 



dq 



■(1 
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)Sp{q) . 



(22) 



ge(g)2 

Note that Eqs. [TO] and [^H involving the static structure 
factor Sp{q), should in principle involve an integral over 
the full dynamical response function Sp{q,uj), see for ex- 
ample ref. [5|. In addition, Eq. l22l omits a term involving 
the square of the excitation energy w, [5|. This term is 
small in the classical limit lu ^ T. We have calculated 
the dynamical response function for neutrons 5'„(g,a;), 
using a semiclassical approximation, in ref. [ll|. How- 
ever, we have not yet calculated the proton dynamical 
response Sp{q,uj). This should be calculated in future 
work and Eq. [22] verified. 

Table [Till presents results for AJL and k. We find that 
K increases some what with density as Z* decreases. To 
gain qualitative insight, we approximate Sp{q) as a single 
peak and write. 



irkpk'^T 
12a2 ' 



■AqZ* 



(23) 



The relatively high thermal conductivity of nuclear pasta 
may keep the pasta phases at almost the same tempera- 
ture as the core of the star. Note that this conductivity 
is from the electrons. A strong magnetic field can reduce 
the electron contribution to the thermal conductivity in 
directions perpendicular to the field. In future work we 
plan to study the nucleon contributions to k. These may 
be important in the presence of a strong magnetic field. 



lar dynamics (MD) simulations of a semiclassical model. 
Our model includes Coulomb repulsion and reproduces 
nuclear saturation. Furthermore, it is directly applicable 
at low densities in a phase with spherical nuclei as well 
as for complex pasta phases. Finally, the model makes 
no assumptions about the pasta shapes. Instead the MD 
simulations could produce any shapes. 

We assume the shear viscosity is dominated by the mo- 
mentum carried by electrons and that this is limited by 
electron-proton scattering. We calculate the static struc- 
ture factor Sp{q) of the protons from our MD trajectories 
for simulations involving 40,000 to 100,000 nucleons. We 
find a peak in Sp (q) that can be characterized by a height 
Z*, a width Ag and a position q* . The height Z* rep- 
resents the effective number of protons that an electron 
can scatter from. This is reduced from the total number 
of protons in a nucleus or pasta cluster because of form 
factor and ion screening effects. 

The shear viscosity can be approximated 



TTTle , q 



20a^'-AqZ* 



(24) 



with Ue the electron density. We find that rj, for nuclear 
pasta, is somewhat increased over that for spherical nu- 
clei, because Z* is much less than the total number of 
protons in a piece of pasta. However rj is not increased 
by orders of magnitude just because of the non-spherical 
pasta shapes. This is in contrast to conventional complex 
fluids where large non-spherical molecules can dramati- 
cally increase 77. For nuclear pasta, the non-spherical 
shapes only impact rj through their influence on Sp{q) 
and the parameters Z*, Aq, and q*. 

Different effective interactions may lead to different 
sizes and shapes for the nuclear pasta. Indeed it may be 
very difficult to predict these accurately. However one 
can calculate transport properties such as rj, the thermal 
conductivity k and the electrical conductivity in terms 
of Sp{q) and the parameters Z*, Aq, and q* . In future 
work we hope to develop a simple way to estimate these 
parameters, given a pasta configuration. In addition we 
plan to calculate other mechanical properties of nuclear 
pasta such as the shear modulus, and extend our MD 
pasta simulations to additional values of proton fraction, 
temperature and density. 
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